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Nuclear model effects in neutrino-nucleus quasielastic scattering are 
studied within the distorted wave impulse approximation, using a relativis- 
tic shell model to describe the nucleus, and comparing it with the rela- 
tivistic Fermi gas. Both charged-current and neutral-current processes are 
considered and, for the neutral-current case, the uncertainties that nuclear 
effects may introduce in measurements of the axial strange form-factor of 
the nucleon are investigated. 
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1. Introduction 

The interest in neutrino-nucleus scattering physics has rapidly increased 
in the past few years, triggered by the rich experimental program aimed at 
studying neutrino properties which is currently under development. The 
proper interpretation of future new data requires an accurate treatment of 
neutrino- nucleus interactions, in order to minimize systematic uncertainties 
due to nuclear modeling. Additionally, the new neutrino facilities provide 
an opportunity to measure nucleon properties, such as the nucleon weak 
form- factors. Such measurements will be obtained from neutrino- nucleus 
scattering data and, again, an accurate treatment of nuclear model effects 
will be needed for a proper interpretation of the results in terms of single 
nucleons. Initially a Fermi gas description of the target nucleus was consid- 
ered to be appropriate, but several studies have shown that this description 
needs to be improved. To this purpose an important guide is provided 
by electron-nucleus scattering, for which both theory and experiment have 
reached a very high degree of sophistication. 

In this contribution we address the problem of nuclear model effects in 
neutrino-nucleus quasielastic (QE) scattering, which has been studied by 
many authors in the past few years ^ |21 El IH El HU • We compare relativistic 
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Fermi gas (RFG) calculations of cross sections with results obtained within 
a relativistic shell model (RSM), originally developed for describing (e, e'N) 
reactions, considering both charged-current (CC) and neutral current (NC) 
processes. The former are relevant for the physics of neutrino detectors, 
while the latter can be used for measuring the nucleon axial strange form- 
factor. In this case, we also consider the ratio of NC over CC cross sections. 



2. Formalism 

Let us consider the QE processes 

Up + A^iT+p+iA-l) (1) 
Vn + A—>Vn+p+(A — i), (2) 

where a neutrino of four momentum K = (e, k) interacts with a nucleus 
A, producing a final state in which a lepton of momentum K' = (e',k'), 
an emitted nucleon of momentum Pjy = (En,Pn) and the (unobserved) 
residual nucleus (A — 1) are present. Following standard procedures the 
exclusive cross section for these processes can be written as a contraction 
of a leptonic and a hadronic tensor 

,,VT ocn^W". (3) 
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The exclusive cross section is then integrated over the momentum of the 
final lepton and/or nucleon in order to obtain the observables of interest. 
In the following, for the CC processes Q we integrate over the emitted 
proton and consider the inclusive cross section (da/dT^), where is the 
outgoing lepton kinetic energy. For the NC processes (J2J) we integrate over 
the undetectable outgoing neutrino and consider (da/dTpj), T/v being the 
kinetic energy of the emitted nucleon. The leptonic tensor in Eq. Q is 
given by 

^ = KpK' v - 9llv K ■ K' + K'^K V - ie^ pa KPK' a . (4) 

The hadronic tensor W^ v is given in general as a bilinear combination of 
matrix elements of the full nuclear weak current, taken between the target 
nucleus ground state and a final state written as a product of the residual 
nucleus (A-l) times the outgoing nucleon scattering state 4>n, and summed 
over all the states of the residual system: 

W^= £ (A-l,^ A r|>(q)|A)(A-l ) ^|i iy (q)|A)*5(^+a;-£; A _ 1 -£; J v). 
(A-l) 
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We calculate W^ v within the impulse approximation, assuming (i) that 
the incident neutrino interacts with only one nucleon which is then emitted, 
while the remaining (A-l) nucleons in the target are spectators, (ii) that the 
nuclear current is the sum of single nucleon currents, and (iii) that the tar- 
get and residual nuclei can be adequately described within an independent 
particle model. Under these assumptions the matrix elements contributing 
to W^ v are greatly simplified and reduce to single nucleon matrix elements 

(A-l^NlPW^ifalJ^B). (6) 

Here Jg N is the single nucleon current operator, which we parametrize in a 
standard way in term of vector and axial weak form-factors (see for exam- 
ple [7j) and ipB and 4>n are the wave functions describing the initial bound 
nucleon and the outgoing nucleon, respectively. In the results presented in 
next section, these wave functions are calculated using the "Madrid-Seville" 
model 8 , originally developed for describing exclusive electron scattering 
reactions and later employed extensively also for studying neutrino scatter- 
ing puna he]. 

In this model the bound nucleon wave functions are obtained as solutions 
of a Dirac equation derived within a relativistic mean field approximation 
from a Lagrangian containing a, u and p mesons. Several descriptions are 
possible for the outgoing nucleon wave functions. In the simplest approach, 
final state interactions (FSI) effects are neglected and 4>N is given by a simple 
plane wave Dirac spinor (plane wave impulse approximation, PWIA). In a 
more realistic approach <f>N is obtained as the solution of a Dirac equation 
containing a phenomenological relativistic optical potential (ROP), obtained 
from fits of elastic proton-nucleus scattering data. Such potential has a real 
part, describing the rescattering of the emitted nucleon, and an imaginary 
part, taking into account the possibility that it is absorbed into unobserved 
inelastic channels. This description (referred to as ROP in the following) is 
appropriate for the calculation of cross sections (da / dT]\[) , where the final 
nucleon is assumed to be detected. On the other hand, when considering 
inclusive cross sections (da/dT^) a selection of the single- nucleon knockout 
channel cannot be made, and the contribution from the inelastic channels 
should be retained [T][I|. Within our approach, a simple way to do so is to 
consider the outgoing nucleon wave functions obtained by setting to zero the 
imaginary part of the ROP (real-ROP approach), thus taking into account 
the conservation of the incident flux. Another possibility is to consider for 
<pN the solutions in the continuum of the same relativistic mean field (RMF) 
equation used to obtain the nucleon bound states. Recent studies of scaling 
properties of inclusive CC neutrino-nucleus QE scattering seem to favor the 
RMF approach [3j. 
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Fig. 1. Inclusive cross sections da/dT^ for the CC processes Q on le O versus the 
outgoing muon kinetic energy T M for E v = 1 GeV. The left panel shows results 
in PWIA and the right panel shows the effects of FSI. RFG cross sections are 
calculated for Fermi momentum kp = 216 MeV. 

We then compare our results with those obtained within the RFG model, 
(see for example [2] ) , including in the latter also a phenomenological energy 
shift (uJsh), which is introduced in studies of inclusive electron scattering in 
order to get the correct position of the QE peak. In light of current /possible 
experiments, we focus here on incident energy E u = 1 GeV and on oxygen 
and carbon targets. 



The differential CC cross sections (da/dT^) for the QE scattering of 
muon neutrinos on ie O are displayed in Fig. ^ as a function of the outgoing 
muon kinetic energy [Q. In the left panel FSI effects are neglected and 
the RFG results are compared with the RSM-PWIA model. We observe 
that, when the RFG energy shift uj s h is taken into account, the differences 
between RFG and RSM-PWIA are very small. The right panel of the figure 
illustrates the effects of FSI within the RSM, by comparing the RSM-PWIA 
curve (solid line, the same as in the left panel) with the curves obtained us- 
ing the real-ROP and RMF approaches outlined above. The sharp peak 
structure observed at large (small energy transfer) is typical when real 
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Fig. 2. Exclusive cross sections da/dT^ for the NC processes (left panel) and 
for the CC processes (Q (right panel). The target nucleus is 12 C and the neutrino 
energy is E v = 1 GeV. RFG cross sections are calculated for Fermi momentum 
k F = 225 MeV. 

potentials are used in the description of the outgoing nucleon. The real-ROP 
and RMF curves, very close to each other, show that FSI produce a reduc- 
tion of the cross section of about 10-15%, which is not too large but may 
have non-negligible effects in the data analysis of experiments measuring 
neutrino properties. 

Let us now consider the differential cross section (da/dT^) for z/„ in- 
duced proton knockout Besides NC processes (left panel), we also con- 
sider CC cross sections (right panel), which will be later used to construct 
ratios of cross sections. We see that, when FSI are neglected, the RSM and 
RFG give results that almost coincide. The inclusion of FSI, which in this 
case is treated using the full complex ROP, produces an important reduc- 
tion of the cross section (~ 50%). The effects are very similar for NC and 
CC processes. 

It is clear that such big effects may have a very strong impact on the use 
of separate NC cross sections for measurements of the axial strange form- 
factor G S A . This is illustrated in the left panel of Fig. |31 where the NC cross 
section integrated over T/v is plotted as a function of g s A = G S A (Q 2 = 0). 
Here a simple dipole parametrization has been assumed for G S A (Q 2 ), with 
the same cut-off mass used for the non-strange axial form-factor (Ma = 
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Fig. 3. Left panel: cross section for NC neutrino induced proton emission integrated 
over the outgoing proton kinetic energy, versus the Q 2 = value of the axial strange 
from factor. The target nucleus is C and the incident neutrino energy is E v = 1 
GeV. Right panel: ratio of NC and CC integrated cross sections (note the smaller 
range of values for g s A ). 



1.032 GeV). 

Is is well known that, in order to extract the nucleon strange form- factors 
from measurements of neutrino-nucleus cross sections, nuclear model effects 
can be largely canceled by considering appropriate ratios of cross sections. 
Several observables have been considered J7JH5], but the most realistic, from 
the experimental point of view, seems to be the so called NC over CC ratio, 
obtained by dividing exclusive NC cross sections by the corresponding CC 
ones. The large cancellation of nuclear model uncertainties in this ratio is 
illustrated in the right panel of Fig- El where the NC/CC ratio of integrated 
cross sections is plotted versus g s A . We see that the very different curves 
shown in the left panel here almost collapse on a single line, making it 
possible to extract the value of g A , provided the experimental errors are 
sufficiently small. 

It is interesting to compare the impact of nuclear model effects on the 
determination of the axial strange form- factor G\, from a measurement 
of the NC/CC ratio, with other possible uncertainties, in particular those 
due to the other form-factors of the nucleon, namely the non-strange axial 
form-factor Ga (assumed to be measured independently in CC processes) 
and the vector strange form-factors (assumed to be determined from parity- 
violating electron scattering data ,7J. This is done in Fig. where the 
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Fig. 4. NC over CC ratio of cross sections dajdT^ for scattering on 12 C at 
E v = 1 GeV, for two values of g s A as indicated. The left panel shows the effects 
of using different nuclear models. The middle panel illustrates the uncertainties 
related to the value of the axial cutoff mass. The right panel gives an example 
of the possible uncertainties due to the magnetic strange form- factor G S M , here 
assumed to have a dipole Q 2 dependence (see text). In the middle and right panels 
we show only RFG results, which are calculated for hp — 225 MeV and uj s h = 0. 

ratio of the NC and CC differential cross sections of Fig. [2] is plotted as a 
function of the emitted nucleon energy T\r. The left panel shows the size of 
nuclear model uncertainties, which are very small. Again, the axial form- 
factors Ga and G S A are parametrized assuming a dipole dependence with the 
same cutoff mass Ma- The middle panel illustrates the effects of changing 
Ma in the range 1.00-1.06 GeV. Finally, the right panel shows the possible 
uncertainties due to the magnetic strange form- factor G S M , for which we 
assume a dipole parametrization Gfj = ^i s /(l + Q 2 /My) 2 with M v = 0.71 
GeV 2 and /i s = 0.37 =L 0.30 lOj, while we assume the electric strange form- 
factor G S E to be zero. At present little is known about the Q 2 dependence of 
the vector strange form-factors, and the results shown here are meant to be 
only an illustration of the impact that the current large errors on them may 
have on the NC/CC ratio. We see that the effects shown in the middle and 
right panels of Fig. 0] are much larger than those due to nuclear modeling. 

In conclusion, we have calculated CC and NC neutrino-nucleus scatter- 
ing cross sections in a relativistic shell model approach, with FSI effects 
taken into account within the relativistic impulse approximation. Our re- 
sults for inclusive CC QE scattering show that FSI effects, although not 
being extremely large, can still be sizable (~ 10%) at the relatively high in- 
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cident energy E v = 1 GeV. For exclusive NC QE processes FSI effects turn 
out to be very large (~ 50%) even at high energy and may prevent a precise 
extraction of the axial strange form-factor of the nucleon form separate cross 
sections. However, these effects are almost canceled when the NC/CC ratio 
is considered. Within our model, nuclear model effects on ratios turn out 
to be much smaller than the uncertainties introduced by the single nucleon 
form-factors and can thus be considered well under control. Of course other 
uncertainties, such as those coming from non-QE contributions to the cross 
sections, or due to other FSI effects not taken into account in the present 
approach, should be carefully considered. 

The results presented in this contribution have been obtained in several 
fruitful collaborations with W.M. Alberico, J. A. Caballero, M.C. Martinez 
and J.M. Udias. 



REFERENCES 



[1] C. Maieron, M. C. Martinez, J. A. Caballero and J. M. Udias, Phys. Rev. C 
68 (2003) 048501. 

[2] W. M. Alberico et al, Nucl. Phys. A 623 (1997) 471. 

[3] W. M. Alberico et al, Nucl. Phys. A 651 (1999) 277; Phys. Lett. B 438 (1998) 
9. 

[4] J. A. Caballero et al, Phys. Rev. Lett. 95 (2005) 252502; J. A. Caballero, 
|arXiv:nucl- th/0604020| 

[5] M. C. Martinez et al, Phys. Rev. C 73 (2006) 024607. 

[6] see for example A. Meucci, C. Giusti and F. D. Pacati, Nucl. Phys. A 
744 (2004) 307; Nucl. Phys. A 739 (2004) 277; B. I. S. van der Ventel 
and J. Piekarewicz, Phys. Rev. C 69 (2004) 035501; C. Bleve et al, As- 
tropart. Phys. 16 (2001) 145; O. Benhar et al, Phys. Rev. D 72 (2005) 
053005; J. E. Amaro et al, Phys. Rev. C 71 (2005) 015501; Phys. Rev. 
C 73 (2006) 035503; J. Nieves, M. Valverde and M. J. Vicente Vacas, 
Phys. Rev. C 73 (2006) 025504; J. Nieves, J. E. Amaro and M. Valverde, 
Phys. Rev. C 70 (2004) 055503; T. Leitner, L. Alvarez-Ruso and U. Mosel, 
|arXiv:nuc l-th/060 1103| 

[7] W. M. Alberico, S. M. Bilenky and C. Maieron, Phys. Rept. 358 (2002) 227. 

[8] J. M. Udias et al, Phys. Rev. C 48 (1993) 2731; Phys. Rev. C 64, 024614 
(2001) and references therein. 

[9] W. M. Alberico, S. M. Bilenky, C. Giunti and C. Maieron, Z. Phys. C 70 
(1996) 463; W. M. Alberico and C. Maieron, |a7Xiv:hep-p h/0210uT7| 
[10] D. T. Spayde et al. [SAMPLE Collaboration], Phys. Lett. B 583 (2004) 79. 



